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Abstract 

For the Burgers equation driven by thermal noise leading asymp- 
totics of pair and high-order correlators of the velocity field are found 
for finite times and large distances. It is shown that the intermittency 
takes place: some correlators are much larger than their reducible 
parts. (Pis'ma v ZhETF, v. 71, n.l, January 2000, in press, English 
Translation in JETP Letters, v. 71, January 2000). 

Intermittency implies strong non-Gaussianity of statistics of fluctuating 
fields. This phenomenon is shown by hydrodynamical systems in a state 
of developed turbulence |jl|, ^. In such far from equilibrium situations 
intermittency appears as prevalence of irreducible parts of some fourth order 
simultaneous correlators over reducible ones. 

As for thermal equilibrium, irreducible parts of simultaneous correlators 
of local fluctuating fields turn out to be of the same order as their Gaussian 
parts even in critical region. This property is inherent in the renormaliza- 
tion group method which takes care of interaction of fluctuations through 
renormalization of local field and effective Hamiltonian 

In the recent paper 0] V.V.Lebedev disclosed that the picture can change 
drastically when we pass to time-dependent correlations of thermally fluc- 
tuating quantities. He found that in the low-temperature phase of two- 
dimensional systems of the Berezinskii-Kosterlitz-Thouless kind different- 
time correlation functions of the vortex charge density may greately exeed 
their own Gaussian part. In the same paper the physical cause of such 
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the behaviour is pointed out: at low temperatures the non- simultaneous cor- 
relation function of every order in vicinity of a given poin are defined by a 
single rare fluctuation. One can conclude from this interpretation that the 
intermittency effects may emerge in equilibrium dynamics of a wide range of 
systems. 

In the present paper I consider one-dimensional velocity field evolving 
according to the Burgers equation with the thermal noise term: 

Ut + uu^ - uu^^ = ^{t,x). (1) 

Here u is the dissipation constant and ^{t,x) is random noise with Gaussian 
statistics and the pair correlator: 

xi)) = -iy(3-'S"{x - xi)S{t - h). (2) 

We will consider u as being vanishingly small. The parameter f3 plays the role 
of inverse temperature, so the simultaneous stationary distribution function 
V[u] has the Gibbs form: 

V[u]=N'exp{~pj^[u]}, J^[u] = j dxu^{x). (3) 

Here M is the normalization constant. The expression: 

{u{t, x)u{t, x')) = {2(3)-^6{x - x). (4) 

following from @ correspond to a total absence of velocity correlation in 
spatially distant points in a given time moment. In the present paper some 
asympotics of various non-simultaneous correlators of the field u{t,x) are 
found. The results obtained here show presence of intermittency effects in 
the equilibrium dynamics of the system (p. 

Dynamical scaling exponent z = 3/2 for the problem (|l|)-(@) was discov- 
ered in the paper using dimensional analysis and utilizing the Galilean 
invariance. In the paper the absence of logarithmic divergencies for the 
spectrum u oc was checked in every order of renormalized perturbation 
theory. Thus the function F2{T,x) =< u(T,x)u{0,0) > has /Sx^/T^ as a 
dimensionless argument. First we find the main (exponential) part of the 
asymptotics of the function F2 at (3x^ /T"^ 3> 1 and z/ — >■ 0. The latter limit 
means that the diffusion cannot set up correlation of velocity in the points 
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and X in a time T. The role of the noise in dynamics on the time inter- 
val (0, T) is also negligible. Thus we can consider u{0,y) as a functional 
of u{T, x) and vice versa. The velocity statistics at the time moment T is 
Gaussian what allows us to represent the non-simultaleous correlator F2 in 
the form: 

The variational derivative Q{t,y) = 6u{t,y)/6u{T,x) for u ^ satisfies the 
continuity equation: 

+ UQy + UyO = 0, (6) 

and the condition 9(T, y) = 6{x — y). This Cauchu problem is solved by the 
characteristic method and we arrive to the expression for F2{T, x): 

F2iT,x) = m-'m,o)) = m-HHx~yiT))(^y%^] ), (?) 



C=0/ 



(see §). Here y{T,() is coordinate of the Lagrange particle started at the 
instant t = from the point (: 

y = u{t,y), 1/(0, C) = C, (8) 

and y{T) = y(T, 0). If u(t, y) is discontinuous, then the equation (|) requires 
a regularization. We use physically evident condition that the particle on a 
shock wave front moves with the velocity of this front. Its formal justification 
starting from a finite small viscosity can be found in |Q. 

The expression (^ tells us that the correlator F2 in the limit being con- 
sidered is determined by a most probable initial velocity fluctuation Mo(?/) 
which, evolving, carries out the particle from the point to the point x in 
the time T. The probabilities of initial configurations are defined by the 
functional (^. The desired optimal fluctuation uq{ii) minimizes T[uq] with 
the condition y{T) = x. Let us show that it is the linear profile: 

uoiy) = ul = x/T - y/T, <y <x, uo{y) = 0, y < 0, y > x. (9) 

First, it is evident that the function Uo(y) must have a maximum at y = 0. 
It is also easy to understand the equality Uo{y) to for ?/ < and y > x. 
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Indeed, difference Uo{y) from zero outside the interval (0,x) does not affect 
the trajectory y(t) but J^[uo] increases. The left edge of the distribution 
u(t, x) for 1/ — >■ will be straight line with the slope a = 1/t. Such the time 
dependence can be checked by direct substitution into Burgers equation; see 
also fl^. At t = T the coordinate of the most rapid particle will be equal to 



X. Coordinates of the other particles from the interval (0, x) will be precisely 
the same. Thus, for the class of initial distributions Uo{y) described above 
the plot of the final function u{T,y) has a form of triangle: 

u{T, y) = y/T, < y < x, u^{y) = 0, y < 0, y > x. (10) 

Now let us note that from the Burgers equation it follows: 

dT[u{t,y)]/dt= -2v j dyul<0, (11) 

what means that: 



J'[uo{y)]>n<T,y)]. (12) 

This inequality becomes strict one even for z/ ^ if shock waves were formed 
during the evolution. Consequently, the minimal admissible value of the 
functional J-' is equal to: 

^[n(T,y)]=xV3T2. (13) 

The value of JF on the function uKy) coincides with (0). The exlusion of 
shocks in the time interval (0, T) justified above makes the expression (^) the 
only possible. 

Probability of the initial fluctuation (^) proportional to exp {—PJ^[uo{y)]) 
defines the exponential part of the asymptotics of the pair correlator F2: 

F,{T,x)^exp(^-^y (14) 

It is worth noting that the multiplier {dy{T,Q /dC)(;=o of the 5-function in 
the formula (|^) vanishes on the configuration (||), but it becomes different 
from zero under small variation of Uo{y). In other words, this factor, along 
with unknown pre-exponential factor in the expression (|T^) as a whole, is 
determined by integration over variations Su of the initial velocity field with 
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respect to u^iy). The essential values of 5u are small comparing with UQ^y); 
the parameter of this smallness is (/?x^/T^)~^. However, the integration over 
Su cannot be reduced to the Gaussian one even in the limit px^/T'^ > 1. 
The point is that at — *• the functional J-'[u] is not analytical on the class 
of initial velocity fields u{y) obeying the constraint y{T) = x. The variation 
6J-' turns out to be of the first order in 6u despite that inequality 5jF > 
holds. J-'[u] can be expanded in functional Taylor series in 6u for 6u <^ u/x 
only. The corresponding analysis will be given in another paper and here we 
restrict ourselves to exponential asymptotics. 

Noting that the initial linear profile tranfers all the internal point of 
the interval (0, x) by the time t = T into the point x we conclude that up to 
pre-exponential factor: 



Fn+2 = (u{T,x)f[u{0,yj) u{0,0) \ ~F2(T,x) ~exp 



/3x^ 
'3T2 



(15) 



Here < yi < y2 ■ ■ ■ < yn- It is obvious that the reducible part of this 
correlator is equal to zero. The same initial fluctuation Mq(|/) determines 
leading asymptotics of the correlator $4 = {u{T, x)u{T, x + ai)u{0, a)u{0, 0)) 
at < a < X and < ai <^ a: 

$4 ~ exp ( 1 > ^i,Gauss ~ exp 1 1 ■ (16) 

Here $4,GaMss designates reducible part of $4. To find $4 as a function of the 
parameter a it is necessary to analyse the evolution of the perturbed linear 
profile. In this case the shock waves formation becomes inevitable what 
makes the problem more difficult. It is worth adding that a-dependence 
of the correlator $4 may be related directly to the probability distribution 
function of velocity field gradients. In |ll|, Q it was shown that the latter 
id defined by forming shocks. Proportionality of asymptotics of high order 
correlation functions to the asymptotics of pair correlator is characteristic 



for turbulent-like problems and in such the context was noted in [|T^. 

Correlation functions of the field u{t, x) may be represented in a form of 
functional integrals (see e.g., fl^). In the present paper such the integrals 
were computed in essence by the saddle-point method with the parameter 
j3x^ /T"^ ^ 1 contained into the object to be averaged, but not inherent to 



the action. This approach goes back to the works of I.M.Lifshits (see in O]) 



Later it was generalized to find high order correlators in equilibrium |16| and 
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strongly non-equilibrium problems |T^, |T^, |T^, ^ ^ ^ . Optimal fluc- 
tuation is called also as an instanton by analogy with quantum field theory. 
In the paper long-time asymptotics of the current autocorrelation func- 
tion has been computed for a disordered contact and the large observation 
time was used saddle-point parameter. 

I am greately indebted to Michael Chertkov and Vladimir Lebedev for 
numerous discussions and advices. I am grateful to Gregory Falkovich and 
Thomas Spencer for backing investigations leading to this work. I would like 
to appreciate Michael Stepanov for useful remarks. 
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